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Abstract 

The theory describing the evolution of perfect fluids coupled to 
a Yang-Mills field in the presence of infinite electrical conductivity 
(ideal Yang-Mills fluids) is proved to be hyperbolic. This result is 
obtained by embedding the theory in a symmetric hyperbolic system 
of equations. 

A divergence formulation of the theory of ideal Yang-Mills fluids 
forms the basis of our analysis. This divergence form constitutes a 
covariant and constraint-free formulation. The infinitesimally small 
amplitude wave structure is derived, and found to be determined by 
a sixth-order polynomial. Alfven waves {5P = 6r = 0) do not, in 
general, exist, and occur in very special cases only. The theory is 
embedded in a symmetric hyperbolic system using a generalization of 
Friedrichs symmetrization procedure. This establishes at once both 
algebraic hyperbolicity (all wave speeds real and less than unit) and 
dependence on initial data to within finite order. 



*This work has been supported by NSF grant PHY 84-04035. 



1 Introduction 



A theory of Yang-Mills fluids has recently been discussed in the limit of 
infinite conductivity, ideal Yang-Mills fluids (lYMF), by Choquet-Bruhat 
[?, ?]. This describes relativistic fluids with a colored magnetic fields. In 
this paper, we concentrate on the mathematical structure of such Yang- 
Mills fluids. The main purpose is to estabhsh well-posedness of initial value 
problems in this theory. Specifically, it will be established 

Theorem 1 Consider relativistic hydrodynamics in which initial value prob- 
lems are well posed by suitable equation of state. For a fluid with infinite con- 
ductivity, well-posedness is preserved in the presence of a Yang-Mills field. 

It is unknown whether a field description for colored magnetic fields is 
physically appropriate to describe the dynamics of high density matter at 
finite temperature, in view of our uncertainty about screening of colored 
magnetic fields. The well-posedness result Theorem 1 above shows that weak 
screening of colored magnetic fields, though uncertain, is permissable. 

We also derive expressions for the infinitesimally small amplitude wave 
structure. While finite conductivity decouples the hydrodynamical and "elec- 
tromagnetic waves," this is no longer true in the singular limit of infinite 
conductivity. 

The analysis is based on two reformulations of the theory of I YMF. 
We begin by obtaining the theory of ideal Yang-Mills fluids in divergence 
form following [?, ?]. This formulation is fully covariant and constraint- 
free. It contains the Yang-Mills solutions by dynamically conserving the 
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original flux-freezing constraints. This formulation Cauchy-Kowalewski reg- 
ularizes the original formulation of lYMF (c/. [?]) and, therefore, defines the 
characteristic determinant. Continuing, we establish well posedness through 
symmetrization. Priedrichs & Lax [?] and Priedrichs [?] introduced a sym- 
metrization procedure for systems of conservation laws in the presence of a 
"main dependency relation" and a certain convexity property. In this pa- 
per, a generalization of their symmetrization procedure is presented which is 
appropriate for the divergence formulation of ideal Yang-Mills fluids. This 
symmetrization procedure allows for the theory of ideal Yang-Mills fluids to 
be embedded in symmetric hyperbolic form. Given this, the standard reg- 
ularity results for quasi-linear symmetric hyperbolic systems apply (see, in 
particular. Fisher & Marsden [?] and referenced therein), which estabhshes 
dependence on initial data to within finite order. 

The organisation of this paper is as follows. In Section 2, the theory of 
ideal Yang-Mills fluids with infinite conductivity is stated. This theory is 
formulated in covariant, constraint- free divergence form in Section 3. The 
characteristic determinant is derived in Section 4. A symmetrization proce- 
dure in the presence of conserved constraints is presented in Section 5. Using 
this symmetrization procedure, ideal Yang-Mills in divergence form is shown 
to allow for an embedding in a symmetric hyperbolic system (Section 6). 

2 Ideal Yang-Mills fluids 

The theory of Yang-Mills fluids describes the evolution of a compressible fluid 
possessing a magnetic field, ha, with values in a non-abelian Lie algebra. In 
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the presence of infinite conductivity, the magnetic field is governed by the 
Yang-Mills' generalization of Faraday's induction law. This description may 
be given in the presence of a curved space-time background with prescribed 
hyperbolic metric gab (thereby considering Yang-Mills test fluids). We will 
do so with metric signature (— , +, +, +). 

The fluid to be considered is assumed to be inviscid and compressible 
with given equation of state. For example, the fluid may be modelled with 
a polytropic equation of state, relating the isotropic hydrostatic pressure, P, 
to the rest mass density, r, as 



Here, K is the adiabatic constant and 7 is the polytropic index. K is con- 
served along the fluid's world lines, except across shocks. For a discussion 
of thermodynamics of relativistic fluids in some generality, see Lichnerowicz 
[?]. Let the velocity four- vector field represent the tangent to the world 
lines of the fluid elements. Then U = {u"', r, P) decribes the state of the fluid 
elements, and is governed by conservation of energy-momentum and baryon 
number: 



Here, T^^ = rfu"'u^ + Pg""^ is the stress-energy tensor for a perfect fluid with 
the specific enthalpy, /. For example, / := 1 + in case of a polytropic 
equation of state (|l|). Notice that (H) constitutes six equations in the six 
unknowns U. 



P = Kr^. 



(1) 



VaTfiU) = 0, 

Vairu^) = 0, 

Co{U) := ^2 + 1 = 0. 



(2) 
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The theory of Yang-Mills fluids in the presence of infinite conductivity 
describes the dynamic interation between a perfect fluid and a (charged) 
magnetic field given by a one- form, ha, with its values in a non-abelian Lie 
algebra Q of dimension A^. Q possesses totally antisymmetric structure con- 
stants Cx^i, (see, e.g., [?]). Elements from Q will be indexed in Greek. We 
thus consider a problem in 6 + 4A^ variables U = (u", ha, r, P). To be precise, 
the Yang-Mills ^-valued 2-form 

Fab := VaAb - VbAa + [A^, Ab] (3) 

in the presence of infinite conductivity, 

u'^Fab = 0, (4) 

may be expressed in terms of a ^-valued 1-form ha as 

Fab = tabcdh^u'^. (5) 

Here, eabcd is the natural volume element on the space-time manifold. The 
condition u'^hc = further imposes uniqueness on this representation (|^) 
of the Yang-Mills field. The evolution of ha is governed by the so-called 
Yang-Mills equations (see [?, ?] for a derivation) 

f V,{/i[V]} = 0, . . 

■ \ c{U) ■= h'u, = 0. ^ ^ 

Here, the gauge invariant derivative Va := V + [A, ■] is introduced so as to 
ensure gauge invariance in (^ (c/. [?, ?]). The conditions c\{U) := h^Uc = 
constitute constraints, one on each Lie-element h^, and may be considered 
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the flux-freezing constraints in Yang-Mills fluids following the nomenclature 
in ideal magneto-hydrodynamics. They impose the condition that the mag- 
netic field be purely spatial in the fluid's rest frame. 

The magnetic field-fluid coupling is established by nonzero magnetic field 
energy and stresses. In the limit of infinite conductivity, this amounts to 
adding the stress-energy tensor T^, 



to the fluid stress-energy tensor, Tj^{U), mentioned above. Here, summation 
over the Lie-index is implicit. In anticipation of what follows, it becomes 
appropriate to define the Yang-Mills magnetic stress tensor := h"-hi, = 
Exhihxa- Notice that T^{U) is entirely determined by H^, as h"^ = H^. 
The evolution of Yang-Mills fluids in the presence of infinite conductivity as 
a test fluid is now described by [?, ?] 



where T^'^U) = Tf(U) + T^{U) is the total stress -energy tensor. Notice 
that (P) constitutes essentially 6 -|- 4iV equations. Of course, (||) must be 
closed by (D and (|). 

It is imperative to observe that (P) constitutes a partial differential- 
algebraic system of equations. Without the constraints c\{U) = co{U) = 
the remaining partial differential equations in (|]) do not impose uniqueness 
on the evolution of the unknowns U. In other words, (^) is not regular in the 




(7) 



' VaT^' = 0, 

< Va(™«) = 0, 



(8) 



cx{U)=0, 
co(U) = 0, 
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sense of Cauchy-Kowalewski (c/. [?]). On the other hand, systems of purely 
partial-differential equations are the most commonly studied, both analyti- 
cally and numerically. It now becomes of interest to first obtain the theory 
of ideal Yang-Mills fluids in covariant form as a system of precisely this type 
with no constraints. 

In what follows we will refer to the theory described by (H) as the theory 
of ideal Yang-Mills fluids, analoguous to the nomenclature ideal magneto- 
hydrodynamics in the case of a single magnetic field {N = 1) with infinite 
'electrical' conductivity. 

3 Ideal Yang-Mills in divergence form 

Consider the initial value problem for ideal Yang-Mills fluids. Maxwell's 
equations require initial data, U^'^\ to satisfy certain compatibility condi- 
tions. Let S denote an initial, space-like submanifold with (time-like) normal 
one-form, z/^, and write 

V, = -Z/a(z/cV^) + (VE)a (9) 

for Va on S. The initial data U^^^ are then subjected to the conditions {of. 
[?] for N = 1) 

f Z/6(Vs)a{/i[V]} = 0, , . 

I c,{U) = 0. ^'"^ 
This constitutes a covariant formulation of the condition that the initial 
magnetic field be divergence free in the space-like, i^a^orthogonal submani- 
fold S, and that the magnetic field is purely space- like for observers whose 
world-lines have tangents u"". 
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Application of the technique described in [?] to each Kx, X = 1, ■ ■ ■ , N, 
yields [?] 

Theorem 3.1 The equations of ideal Yang- Mills fluids including the flux- 
freezing constraints can be stated as 

VaT'^' = 0, 

Va{hK''^+g-'c{U)} = 0, , . 

Vairu") = 0, ^ ^ 

This system must be closed by 

VaAb - VbAa = eabcdh'u'^ - [Aa, A^]. (12) 

The vector is any prescribed time-like vector field. This system is equiv- 
alent to ^ in regions where the flow is continuously differentiable. The 
standard jump conditions across surfaces of discontinuity for this system are 
those of conservation of energy-momentum, baryon number and Yang-Mills 
equations. 

Proof : Clearly, we need only show that a solution to an initial value 

problem in the new formulation with Cauchy-data satisfying compatibity 

conditions ([T0|) on an initial, space-like hypersurface S, yields a solution to 

the original system of PDAE's (H). We will do so by showing that cx satisfies 

the a homogeneous wave equation with vanishing Cauchy-data: 

' Dc = in 
c = on E, (13) 

v^VaC = on S. 

Here, we have defined the Yang-Mills wave operator □ = V'^Vc+[(V'^Ac), ■]-!- 
[A^Va, ■] + [v4^, [Ac, ■]]. -D^(S) denotes the future domain of dependence 



of S, and z/^ denotes the normal one-form to S. In what follows, eabcd = 
y/—g[abcd] will denote the volume element on our four-dimensional manifold 
with = det[gab]- Here, [abcd\ is the totally antisymmetric symbol such 
that [0123] = 1. 

Ricci's identity implies the identity [?] 

VaVhi2hK'^) + ^[F,,, F,, - t.defh^ufy^'^'' = 0. (14) 

By (|1|), Fab = eabcdh^W^ is enforced, so that V aVb{'2h^''v!''^) = 0. Conse- 
quently, we are left with the Yang-Mills wave equation for c: 

= Va^biuj^') = V,Vfc(/i[Vl) + (V^Ve)c = DC. (15) 

On S, we may write 

V„ = -Z/a(z/^Ve) + (Vs)a. (16) 

Using this, Cauchy-data satisfying (|1^) yield 

=u\W''{hyaU^+ gabC)} 

= -u'u''{u'^V,){h[aUb]} + V^(Vs)n^[a«b]} + ^'^bC (17) 

because h^aUf,] is antisymmetric. 

This forces c = in D^(S), and the proof is complete. □ 
We remark that for isentropic fluids, where the entropy is constant ev- 
erywhere, the fluid variables r and P depend on essentially one parameter, 
in view of dP\s = rdf. In this event, conservation of energy-momentum, 
VaT"'''{U) = 0, together with conservation of baryon number, Va(ru") = 0, 
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may be seen to yield conservation of the constraint cq{U) = 0. (0) then 
reduces to 

V^T-^" = 0, 

Vaj/i'V] + g''''c{U)} = 0, (18) 
Vairu'') = 0. 



(11) strictly contains ideal Yang-Mills solutions to PDAE (M). Ideal Yang- 



Mills solutions are characterized by vanishing of the flux-freezing constraints 

cx{U)=0. (19) 

Initial data satisfying the compatibility conditions ([10|) ensure that the c\{U) 
are preserved as zero (as in the case of MHD [?]). The formulation of the con- 
straint co{U) in terms of the last equation in ( [iTl) is taken from the divergence 
formulation of ideal MHD [?]. 

4 Ideal Yang-Mills' characteristic determinant 

Infinitesimally small amplitude waves may be defined through the normal 
cone, A/", of the one-forms Ua normal to their time-like surfaces of propaga- 
tion. We will consider waves whose velocity of propagation is strictly less 
than unity, i.e., waves with u'^Uc 7^ 0. The condition on Uaef/ is given by 
the vanishing of the characteristic determinant (pointwise on the space-time 
manifold) 

D{U;u,) ■.= det-^ = 0. (20) 

Here, B indexes the unknowns U = [u"" ,ha,P,r), as defined before. Clearly, 
D{U\Va) constitutes a 6 + AN degree polynomial in Va. The wave struc- 
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ture may be classified by finding the factorization of D{U; Va) over miminal 
polynomials in Ua- 

Using Theorem 3.1, we have 

/ K Ri ■■■ Rn \*\ 



det 





A 











A 







(21) 








A 







* 


••• 





* I 





V 

Here, we have introduced the 4 by 4 matrices (when written in component 
form) 

(r/ + /i2)(MV,)5^ + (r/ + /i2)n''z/„ 
:-h>iiya+ihlu,)6', + iy'h^a, 



T3b dT'=''iyc 



(22) 



dhl 



where 



h°'u\ — h^u'^ + g°'''cx{U). The main step towards factorization 



is contained in 

Proposition 4.1 The Yang-Mills characteristic determinant D{U]v a) con- 
tains the factor |y4|^~^. 

The proof of the Proposition is based on a suitable rewriting of the char- 
acteristic matrix [F^]{U\ Va) in ( PtlD and (j2T|), and the observation 



Lemma 4.1 Consider a rank- one update, Mx, to an N x N matrix, M = 
[m^], through an N -dimensional column vector, u^, and row vector, Va, of the 
form M + X[u^Va]- The determinant of Mx = M + X[u^Va] is linear in X. 
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Proof : If either u'' or Va is zero, there is nothing to prove. Let, therefore, 
0. We have 

det Mx = det [[mj + Xu'^Vi], • • • , [m% + Xu'^vn]] 

= det [K - ^^fe], • • • , H-i + Xu'vk-i], [mi + XuW], . . 

K+i + A«V+i],---,[<-t<]] ^ ^ 

= det [K], • • • , [mti], [ml + AuV], K+J, • • • , [m%]]. 

Because of the definition det [H^] - Y.ai-^Y^'' '^'' Hl^Hl^ - ■ ■ H^^ for the 
determinant of an iV x iV matrix [i^^], where e is the iV— dimensional totally 
antisymmetric symbol and the sum is over all iV-permutations cr, it follows 
that each term in our determinant contains A*^ with A; = 0, 1. □ 



Proof of Proposition 4-1-' The result will follow in several steps. The first 
two steps concern rewriting of the characteristic matrix of Yang-Mills' in 
divergence form through column block operations by (unitary matrix) mul- 
tiplications from the right and row block operations by (unitary matrix) 
multiplications from the left. The final result follows from further inspec- 
tion of a resulting 6 by 6 reduced matrix. 

Step (a): Operate on the leading column block through multiplication by 



matrices with unit determinant from the right: 

/id 



D{U-Ua)^[ 
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-A-^S^ id 







id 










id 






(24) 



|id / 



12 



( K- RiA-^S^ R 



R 



N 



\* \ 










A 




|o 


det 






A 




|o 










A 


|o 




I 


* 


••• 





* 



(25) 



Proceeding in this manner to cancel each in the leading column block 
through subsequent multiplication of the matrix in ( p5|) from the right by 
matrices with unit determinant of the form 

/id ■■■ |0 \ 




(26) 



we obtain 



( K- RA-^S Ri 



R 



N 



I* \ 










A 




|o 


D = det 







A 




|o 











A 


|o 




I 


* 


■ ■■ 





* 



= det[ 



(27) 
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Here, we have used 



RA-^S = Y,RxA-'^S\ 



(28) 



Step (b): By a procedure similar to that used in Step (a), but now row- 
wise, we may proceed to cancel each Rx in the leading row block through 
multiplication by matrices with unit determinant from the left. Beginning 
with cancellaetion of Ri, 





( 


-RiA-' 


•• 


■ 


|o \ 







id 






|o 


D{U;Ua) =det 







id 




|o 















id 


|0 




I 










id) 


( K- 


RA-^S 


R2 ■ 


• • Rn 


1* 








A 


|o 


det 





A 














A |0 



V * ••• I* y 

we may continue in this manner until all R\ have been removed: 

/ K-RA-^S ••• I* \ 








A 


|o 


D{U;i^a) = det 





A 














A |0 



I* / 



(29) 



(30) 



(31) 
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It follows that 

( K-RA-^S I* 



D = Ul^det 



ArdetM(f/;z/J. (32) 



V 



Step (c): It may be readily verified that 



A-' = -(mVJ-i[(7^ + - ly-^iy'iya] 
detA = u^iu^Ur?. 



(33) 



Multiplication of [D^] {U ; Va) by an appropriate matrix with unit determinant 
from the right, whereby multipling the leading column block by A from the 
right, leads us to consider 

{KA- RA-^SA I* 



D{U-Ua) =|A|~-Met 



*A 

|Ar-Met[5f](t/;z/J. 



(34) 



Step (d): It remains to ascertain that det[Z)^](f/; Va) constitutes a poly- 
nomial in Va- This is not immediate, in view of the presence of A~^ in [D^ 
and the factor (■u'^z/c)"^ in (^). As it turns out, any rational factor thus in- 
troduced in det[Dj^]{U ; Va) will fortuitously be cancelled. This will be shown 
below. 

Explicit evaluation of A~^S^A gives 

A~^S^A = {u'u,)-^[{u^h^^ - (/i^V,)z/^)mJ + Q\ (35) 

where = {h^'^i'c)id + [u^{{u'^i'c)h^ — {h^'^i'c)ua] represents the remaining 
matrix whose entries form polynomial expressions in Ua. This gives 

RA-^SA = v\u'vc)'^[{h^v^ - h^)ua] + i?Q, (36) 
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where RQ = J2xRxQ^ and h'' = Y.\{hli^c)h^^ ■ This shows that {u^Vc)~^ 
enters in a rank-one update to the [ua- polynomial) KA — RQ, and, therefore, 
introduces a factor {u'^v^^ with k> —I'va each term in detfDf ]([/ ; Va). This 
follows by Lemma 4.1. More elementary, we may introduce a special frame of 
reference in which Ua = (—1, 0, 0, 0). (u'^Uc)^^ then simply enters in the first 
column of [D^]{U;i'a), and the forementioned observation readily follows. 
Inspection of the full matrix [Df]{U; with U = (^^^ P, r), 

f KA- RA-^SA \qp{u'Uc)[u^] + [u^] \qr{u'iyc)[u^] \ 



r[u,Al] |0 l(iiVe) 

V (ei^c)M^] |o |o J 



,(37) 



where qy = dq/dy, q = rf + h?, further shows that the last column is propor- 
tional to {u^Vc)- This dependency introduces a factor {u'^Pc) in det[Z)^](f/ ; z/^) 
which, combined with the foregoing, forces each term in det[D^]([/; !/„) to 
contain {u'^v^)^ with A; > 0: <\ei[D^{U ; v^) forms a polynomial expression in 
Va (homogeneous of degree 10). 

Together, the results from Step (c) and Step (d) show that lAI-^"-^ is a 
factor in the the 6 -|- 4iV-degree homogeneous polynomial D(U ; !/„). □ 

The full factorization is now contained in 

Theorem 4.1 Yang-Mills' characteristic determinant is of the form 

D{U;^a) = -iCi^c){u'^cy''-\i^'i^cfY{U;Ua), (38) 
where F(t/;i/a) is a sixth- order polynomial, homogeneous in Va- 
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Proof : Inspection of [D^]{U; Va) as defined in the proof of Lemma 4.1 
shows that all its entries are even in except for the lower row 

[r[v,A%Q,{u^v,)] (39) 

stemming from conservation of baryon number, V airu"") = 0. det [D^]{U ; z/q) 
is therefore odd in u*. In view of \A\ = (z/^z/c)(m'^z/c)^, this forces D{U; Vg) to 
be odd in Now notice that appears only through {u'^v^ in the scalar 
field D{U; Vg)-, since u^h^ = for all A and u^u^ = —1 (see also Proposition 
5.1 in [?]). The roots of D{U; z/^) are invariant quantities, necessarily invari- 
ant under a sign-change —vf'. It follows that an odd factor of (u'^i'c) 
must factor in D{U; z/^), leaving all other factors even in {u'^v^. The factor 
(u^z/c)^^~^ introduced by therefore, forces {u^UcY^^^ to be a factor 
in det D{U ; z/^). 

The factor {v^v^^ = z/^^ in (|55|), rather than promised by 

|y4|^~^, follows from the fact that the 6 x 4N dimensional elements 

yUfceC, form left null elements of the characteristic matrix [F^]{U : z/^), when- 
ever z/^ = 0. This follows from z/^ being a null element of both A^^ and S'^'', 
i.e., 

vAl = ^cS^' = 0, (41) 

whenever z/^ = 0. This completes the proof. □ 

-'^In [?] the characteristic determinant of ideal MHD [N = 1) is even in , as we work 
there with u^V aS = (instead of V a{{u^ + 1)^} = 0). 
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We conclude that the null cone M at each point on the space-time mani- 
fold is described by 

(1) Entropy waves: u'^u^ = 0, 

(2) Yang-Mills waves: Y{U; Ua) = 0. ^ ' 

Thus, the Yang-Mills waves are completely defined by a sixth-order polyno- 
mial Y{U]Va). 

Recall that in the theory of ideal magneto-hydrodynamics (MHD) we 
find the hydrodynamical waves and the Alfven waves through the zeros of a 
sixth order polynomial, which factorizes over a fourth-order (hydrodynamical 
waves) and a second-order (Alfven waves) polynomial [?, ?]. However, one 
may not expect such factorization to persist in the more general case of a 
Yang-Mills fluid. 

Expressions (27)- (28) show that the Yang- Mills magnetic stress-tensor 
if^, as defined above, represents the entire magnetic field haeQ in the in- 
finitesimal wave structure. It follows that the Yang-Mills waves Y{U ; ^a) = 
are completely determined by {u"' , H^^, P, r) . Because H^b is real-symmetric, 
the Spectral Decomposition Theorem allows for the representation 

Ha, = i:s^sP, (43) 

i=l 

the s^*-*, i = 1,2,3, being three one-forms on the space-time manifold. Since 
Habu"" = 0, it follows that N — 3 without loss of generality in the structure 
of the infinitesimally small amplitude waves in ideal Yang-mills fiuids. 
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The new wave features in Yang-Mills fluids in ideal Yang- Mills fluids are 
contained in the normal cone as defined by the zeros of the Yang-Mills char- 
acteristic polynomial Y{U ; Ua). Of foremost importance is to establish hyper- 
bohcity, i.e., to estabhsh that Y{U; Va) = contains a full set of three (real- 
valued) normal conical sheats which are everywhere space-like. Of course, we 
need not require strict, algebraic hyperbolicity, that is, some of these normal 
sheats may (partially) coincide, associated with multiplicities of the roots of 
Y{U] Vfj) — 0. This may be studied through symmetrization. 

5 On hyperbolicity of systems of conserva- 
tion laws 

Friedrichs & Lax [?] and Friedrichs [?] conceived the deep and general result 
that a system of conservation laws is well-posed when the system (implic- 
itly) contains a conservation law for a convex quantity. A quantity of this 
kind naturally enforces an integral bound on the solution for all time. More 
specific regularity results for finite time are obtained following symmetriza- 
tion, as obtained by such conserved convex quantity. The resulting system 
is in symmetric hyperbolic form ([?]), which establishes at one stroke both 
hyperbohcity and well-posedness. In our study of the Yang-Mills equations, 
a slightly weaker version of Friedrichs' symmetrization procedure is formu- 
lated in view of the conserved, algebraic constraints in the theory of ideal 
Yang-Mills fluids. 
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5.1 Priedrichs' symmetrization procedure 

Ideal Yang-Mills as given in Theorem 3.1 is a particular case of a system of 
M conservation laws of the form 



<aB 




(44) 



Friedrichs [?] introduced certain Properties CI and CII, discussed below, 
which are sufficient for (^4]) to constitute a symmetric hyperbolic system, 



A'^^^{V)VaVA = A'^''iV)VtV + A"^^(V^)V«\4 = f'iV), (45) 



characterized by symmetry of the A"-^^ (in A and B) and positive definiteness 
of A^"^^. Systems of this type are known to give rise to well-posed initial 
value with sensitivity on initial data of finite order. In the context of general 
relativity, the requirement that A"-^^C,a is nonsingular for all time-like one- 
forms C,a furthermore ensures that no infinitesimally small amplitude wave 
exceeds unit velocity. 

Friedrichs' [?] symmetrization procedure applies to systems possessing 
what Friedrichs calls a "main dependency" relation of the form 



Here, Wa is a nontrivial vector and 6 denotes a total variation (this property 
is called CI' in [Friedrichs, 1974]). Evidently, this requires Wa/^ = 0. Then 
differentiation with respect to any yields 



CI: WaSF^^ = 0. 



(46) 



OWa dF"^ 
dvc dV^ 



V.r^ + Wa 



D 



0. 



(47) 
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This shows that for each a 

dWA dF"^ 



(48) 



is symmetric (in C and D). If, furthermore, is such that for some time-hke 

CII: SWA5F^^^^ = 5V^i^^)6V^ (49) 

is positive definite (for all SV^) , then (^11) naturally constitutes a symmetric 
hyperbolic system in Va in the form of 

dWAdF^^ r,_dWA.A (.r,. 

Writing = W aF"^ - it follows that solutions to (|0D satisfy 



WaZ^ = 0, ^ ^ 



where the second equation is a consequence of Property CI. If Property CII 
is satisfied as well, {Wa, can be readily seen to be nonsingular, so that 
any solution to (^), is, in fact, a solution to (44). 



This summarizes the symmetrization procedure constructed by Friedrichs 
[?] for systems possessing a single dependency relation (^6]). It should be 
mentioned that Friedrichs [?] also treats the case of more than one depen- 
dency relations. For the purpose of symmetrization of ideal Yang-Mills in 
divergence form there exists precisely one, main dependency relation, al- 
though in a slightly weaker form than CI as stated above. Not suprisingly, 
this is due to the fact that the divergence formulation allows for a larger class 
of solutions, which includes the nonphysical solutions with ca 7^ 0. 
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5.2 Symmetrization in the presence of conserved con- 
straints 

Ideal Yang-Mills in divergence form contains the constraints ca = as con- 
served quantities when initial data satisfy certain compatibility conditions. 
The system itself allows for solutions with ca 7^ as well, in response to 
nonphysical initial data. It is therefore not suprising that in this generality 
ideal Yang-Mills in divergence form does not satisfy conditions CI or CII. 
Rather, we require a main dependency relation of the form 

YI: Wa5F''^ = 5z'' (52) 

for some nontrivial Wa and some vector field z'^. We will later require that 
z"" vanishes for solutions which satisfy the constraints. Clearly, Property YI 
is weaker than Friedrichs' Property CI. However, the identity 

instead of (^Tf), clearly does not change the symmetry of the coefficient ma- 
trices (^8|). In what follows, nonsingularity of {Wa, will be assumed, so 
that (|5TD enforces equivalence of (^) with (^^. 

For solutions with ca 7^ 0, different from the Yang-Mills solutions ca = 0, 
we may not expect convexity property CII to hold. We do, however, wish 
solutions to ideal Yang-Mills satisfying the constraints 

Co = 0, (54) 

CA = 0, (55) 

to satisfy a symmetric hyperbolic system. To this end, a slightly weaker 
formulation of CII will be used. 
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Definition 5.1 Consider a scalar Q = Q{U^ , 6U^') and a set of scalar con- 
straints, CkiJJ^) = 0. The variation 6U^ 7^ satisfying 5ck{U) = will be 
referred to as the constraint variation of If-"' with respect to the Ck = 0. 
If Q > with respect to the constraint variation of , then Q is said to he 
constraint positive definite with respect to the = 0. 

Notice that for each the constraints in Definition 5.1 define a hnear 
subspace of variations . 

We now define constraint positive definiteness of 

YII: Q:=5WA5F^^ia = 5V^{^'-^)5V^ (56) 

for some time-hke S,a as Property YII. Clearly, YII is weaker than Friedrichs' 
Property CII. 

Property YI and YII are sufficient for to imply a symmetric hyper- 
bolic system in Va. To see this, we use the following construction. 

Lemma 5.1 Given a real-symmetic y4e£(R"',R") which is positive definite 
on a linear subspace V C R", there exists a real-symmetric, positive definite 
v4*e/:(R",R") such that 

A*y = Ay{yeV). (57) 

Proof : Let V"*" denote the orthogonal complement of V. If V"*" = {0}, 
we may take A* = A. Otherwise, pick xeV'^ with ||x|| = 1, and form 

Ai = A + fxxx^, (58) 



23 



where denotes the transpose of x. Ai is manisfestly symmetric. We will 
choose /i so that Ai is positive definite on the space Vi of all vectors of the 
form z = y + Ax, where ytV and A is a real scalar. 
To this end, evaluate 

z'^Az = y^Ay + 2Xx'^Ay + X^x^Ax + A^^ 

> c\\y\\^ - 2|A|M|||/|| + A2(/i2 _ m) (59) 
= (^^-M)(|A|-™)^ + c|MP-^|MP. 

Here, c > is a coercivity constant for A on V, i.e., y'^Ay > c\\y\ \ (yeV), and 

M = \ \A\\. For sufficiently large fx > M'^^'^, we have 

C||y||2__il!^||^||2> ^||^||2. (60) 

Writing A = a||?/||, ( |59D yields 



z^Az>K{a)\\z\\^ (61) 

where K{a) = (1 + a^r'{{fi^ ~ M){a ~ + f }. Clearly, inf,,RK(a) = 

ci > 0. 

We have now proved that there exists an embedding of A |v in a symmetric 
Ai which is positive definite on Vi. 

Repeating the construction above k = dimV^ times, we may exhaust 
V"*", thereby arriving at a symmetric A* = Ak which is positive definite on 
Vfc = R". This completes the proof. □ 



Consider a system (^) which satisfies Property YI, thereby arriving at a 
system (^), which may, alternatively, be written as 

A'^''{V)dtV + A°^^(1^)9,\/a = F{V). (62) 
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Here, now contains the contributions from the connection symbols as 
well. Assuming that (Q) satisfies Property YII with = {dtY, we may 
form a positive definite (A*"^^)* following Lemma 5.1. Indeed, for V such 
that Cx(y) = for all A, Property YII implies A^"^^(y) is positive definite 
on V = PIaI'^^I^^^^c = 0}. According to Lemma 5.1, we may insist 

= A'^^'y (63) 

for all constraint variations yeV. 

We are thus led to consider the symmetric hyperbolic system 

{A'^y{V)dtV + = F{V). (64) 

The arguments above show that solutions to the original divergence sys- 
tem of equations (0) in response to initial data compatible with (p!0| ) (which 



have the property that c = throughout) are solutions to the symmetric 
hyperbolic system of equations ([641). It naturally follows that these solutions 
depend continuously on their data as described in [?]. 

Friedrichs's symmetrizations procedure and the symmetrization proce- 
dure in the presence of conserved constraints as presented in this Section are 
summarized in Table I. 

It would be of interest to also use the symmetric hyperbolic system of 
equations ( |64D for existence proofs of solutions to the divergence formulation 
of ideal Yang- Mills fiuids. To this end, solutions to initial value problems for 
the symmetric hyperbolic system of equations ( |64D in response to data com- 
patible with (p!0| ) must be shown to have the property that c = throughout. 
However, this falls outside the scope of this work. 
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Friedrichs ' procedure 


In this paper, with conserved c = 


Properties Cx 


Properties Yx 


CI: WaSF^"" = 
CII: SWASF^'^^a 

positive definite 
(CIII when more than one 
dependency relation [?]) 


YI: Wa^F"^ = Sz^ 
YII: 6WA6F^^^a 
constraint positive definite 
YIII: c = if (|10D satisfied 


dWA dF'^'^'rj T/D _ dWA fA 


fdWAdF*'^\*fi'[/D 1 dWA dF°'^ ft -[/D SWa fB 



Table 1: Symmetrization procedures: summary of Friedrichs' symmetriza- 
tion procedure and the symmetrization procedure in the presence of con- 
served constraints, as presented in this paper. Friedrichs' procedure requires 
a system VaF""^ = to possess Properties CI and CII (in the presence of a 
single dependency relation), and the procedure presented here requires the 
system to possess Properties YI, YII and YIII. 

6 Ideal Yang-Mills in symmetric hyperbolic 
form 

We will show that ideal Yang-Mills satisfies a symmetric hyperbolic system 
of the form (^5]). Using the symmetrization procedure in the presence of 
constraints as outlined in the previous section, this result follows by showing 
that ideal Yang-Mills in divergence form (|Tl|) satisfies Properties YI and 
YII. The proof of Theorem 3.1 establishes that (0) satisfies Property YIII. 
It is well-known that the equations of relativistic hydrodynamics can be 
written in symmetric hj^erbolic form [?, ?]. The desired result, therefore, 
will be established "bootstrap-wise," by showing that ideal Yang-Mills can be 
written in symmetric hyperbolic form, whenever the equations of relativistic 
hydrodynamics can be written in symmetric hyperbolic form. 
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For the purpose of symmetrization, we use a second divergence form of 
ideal Yang-Mills fluids (following [?]), in which the equations of hydrody- 
namics are more similar to those already discussed in the literature: 



VaT'^' = 0, 

Vairu'') = 0, ^ ^ 

Va(r5u") = 0. 

This system may readily be seen to be equivalent to ( pT] ) for continously dif- 
ferentiable solutions in view of the thermodynamic relation dP = rdf—rTdS. 
Note, however, that the jump conditions across surfaces of discontinuity as 
follow from a weak formulation for this system (|65|) are improper (see [?] for 
a discussion on this point in magneto- hydrodynamics). We organise (|65D as 
follows 



where 



VaFf + VaF:f = ff + (66) 
V„Tf, 

VaFf = \ V„(™-), // = 0, 

V.{e(«' + 1)}, (g7) 

( TpaA\ — J ^a-^m ) {B _ j 

^ay-l^m A — 1 Tj , ,ab Jm ) 



m ) 



Here, ca^j/ are the structure constants of the Yang-Mills Lie algebra, Q. As 
mentioned before, these equations must be closed by ([T2|) , 

VaAk - VbAa = Fab ' A] = Gab- (68) 
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Before turning to the symmetrization of the nonhnear system of equations 
(|65[), we wish to remark that with temporal gauge the closure equations can 
be seen explicitly to pose no difficulties in the symmetrization process. To 
see this, let is any given time-like vector field, and consider the equations 

(rVc)Aa - rVa^c = rCca. (69) 

Because Vacu"'' = enforces V[a-Ffec] = iaFbc] + [A[a, Fbc]] = 0, and [Fiab, A^] = 
V[a[Ab, Ac]], the right hand-side in (|68D satisfies V[aGbc] = (using the re- 
sult that c{U) = 0). Denoting with Greek subscripts contractions with 
orthogonal vectors, (^) thus may be seen to ensure conservation of 
VaAjS — VpAa = Gap- Eqs. (11)-(12) therefore impose (75) as evolution 
equations for Aa- With conservation ^^Ac = 0, any solution to (p9| ) satisfies 

{edc)Aa = CGca + AcdaC (70) 

with coefficient matrices 

(eVc)^?^ (71) 

For this reason, the closure equations (0) with temporal gauge need no fur- 
ther consideration in the symmetrization of the equations of ideal Yang-Mills 
in divergence form. 

In the variables (m*, h^, f, S) (^) is, of course, regular in the sense of Cauchy- 
Kowalewski (by its equivalence to (pi])), and possesses therefore no depen- 
dency relation. If we consider (p5| ) in = (f^, h'', /, S), however, where Va is 
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a three-parametrization of such that u"^ = —1, this system becomes overde- 
termined by one and a main dependency relation results (as in Property YI). 
This will be made explicit below. 

6.1 Symmetrization of hydrodynamics 

The equations of relativistic hydrodynamics, 

VaFf = 0, (72) 

as follow from ideal Yang- Mills in divergence form with vanishing Yang-Mills 
field, may be written in symmetric hyperbolic form following the arguments 
given by Priedrichs [?], Ruggeri & Strumia [?] and Annile [?]. Priedrichs [?] 
emphasis that symmetric hyperbolic systems will not generally be in covari- 
ant form; the construction of systems of this type is intended as a method 
of establishing well-posedness, an aspect of which is frame independent. Por 
symmetrization of the equations of relativistic hydrodynamics, the veloc- 
ity four-vector, u'', may be represented by a three-parameter representation, 
u'' = u''{va), in view of cq = u'^Uc+1 = 0. Consider = j;{ua, f — TS, T) as 
a function of V/j = {va, T, /). Notice that variations 6V^ satisfy Su^ = by 
construction. In this parametrization, both Priedrichs' Properties CI & CII 
hold true. Indeed, we have a main dependency relation and an associated 
quadratic, Q/, of the form [?] 

WiSFf^ = Ub6Tf + (/ - TS)6{ru'') + T5(r5u") = 0, 

Qf.^SWASFf^ = T{5[f]5Tf + S[^]5{ru'')} + fwi6F^^ (73) 
= T{5[^]5Tf + S[i^]5{ru-)}. 

Ruggeri & Strumia [?] and Anile [?] have shown to be positive definite, 
provided that the following conditions on the free enthalpy, G, and the sound 
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velocity, a, are satisfied [?]: 

f the free entlialpy: - 
I tlie sound velocity: a 




(74) 



We will now proceed with the full equations of ideal Yang- Mills in divergence 
form. 

6.2 Symmetrization with Yang-Mills field 

Ideal Yang-Mills in divergence form will be shown to be satisfying Property 
YI and YII, provided that the hydrodynamical conditions (ff^) have been 
met. 

Property YI. Recall from (5) and (14) the contribution to the stress- 
energy tensor by the magnetic field, T^, and the new Yang-Mills tensor, u;"'', 
respectively, given by 



Let us consider Wa = j^{ua, ha, f) with parametrization Va = {va, ha,T, f). 
The total variation, 6, then satisfies co{U) := + 1 = by construction, as 
before. The variational expressions 



(76) 



(75) 



+g°-^hM'' - h''6h'' - h'^dh''} 

-h^du"" - u'^ihcSh") - h^iucSh^) - cSh", 



(77) 



introduce the identity 



nb 



26{ch^) = Sz". 



(78) 
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This establishes that Property YI is satisfied: 

WA5{Ff + F^^) = bz\ (79) 

Property YII. Ideal Yang-Mills in divergence form can be shown to be 
satisfying Property YII by consideration of the extension to relativistic hy- 
drodynamics as introduced by the Yang-Mills fields. 

Proposition 6.1 // the equations of relativistic hydrodynamics are convex, 
then ideal Yang-Mills in divergence form is constraint convex. 

Proof : Let Qf denote the quadratic variation associated with the equa- 
tions of relativistic hydrodynamics as in (^), and suppose conditions (^) 
have been met. 

Consider the variation, 5, introduced by the total variation of Vc {i.e., 
variations satisfying 5cq{U) = 5{v? + 1) = 0). The variational expressions 
5u,5T^U = h\u%,){6uM'f 

5h,5uj-^ia = m,){5hM') + {ic5h%u,5h-) ^^^^ 
-{iM%K5h'^) - {u%){Sh,5h'^) + {^M')Sc, 

introduce the identities 

Su,6T:i^'^a-Sh,6uj''% = {u%,)[h\6uf + {6hf] + 2[{^M%hJh'^) 

-{h%,){6ujh')] - 2{^Jh')Sc = Qo. ^ ^ 

If 6h^ varies freely, disregarding 6cx = 0, the quadratic Qq fails to be positive 

definite. 

In the presence of the constraint variations with respect to the cx = 0, 
the last term in the expression for Qq in (|8TD vanishes identically: 

- 2{CcSh')Sc = ^ ^2{^Ml)Scx = 0, (82) 

A 
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and we are left with 



Qr 



ab c 



(83) 



This expression is (constraint) positive definite in constraint variations with 
6h 7^ 0. To see this, it is convenient to work in a local Lorenzian frame of 



reference g 



ab 



ab 



dia(-l, 1, 1, 1) with {x"} 



a 



1, 2, 3, such 



that u'' = (1,0,0,0). The magnetic field, h^, and the variations 6u^ and Sh\ 
may then be expressed in terms of three-component, spatial vectors, k^, 6v^ 
and respectively, as 



Su'' 
5h\ 



/3' 



(84) 



Notice that 5e\ follows from 5c\ = 6ex + k'^6v^ = 0. In the constraint 
variation of 6, the right hand-side in the quadratic variation in (|8l|) becomes 
a A— sum of quadratic forms SWaxQx^^Wbx in SWax = {Sva,Skax), where 



Q 

With {u%, 



AB 
A 



(85) 



^0 > 0, the six eigenvalues, s;, of i^ay be evaluated as 



si = 0, 

S2 = {U%c) > 0, 



A±{A'-[{u%,fe-{k%,f]y^>0, 

A±{A^ -k^\e\y^ >o. 



(86) 



where A = ^{u''C,c){k'^ + !)• Here, we have used C,"^ = — + ^"^a < 0, so 
that ^2 = |^2| ^ ^"^^ > |^2| ^2 ^ ^2 ^^^Yi a = 1, 2, 3. Notice that the 
eigenvector of Qf^ associated with si = is (/c^, 0). Therefore, WaxQx^Wbx 
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is positive semi-definite, and positive definite whenever Sk^ ^ 0. The sum 

T.^AxQt'^WBx, (87) 

therefore, naturally shares the same properties. 

li Qf is positive definite with respect to SV^, the preceeding result shows 
that 

SVASU'' = Qf + Qm (88) 

is constraint positive definite with respect to the Ca = 0. □ 

The constraint symmetrization procedure from Section 5.2 may now be 
applied to ideal Yang-Mills in divergence form with Va as given above. It 
remains to ascertain the orthogonality requirement 

= {u'^h^,)cx^,h''-A>^ - {h^h''^)cx^,{A>^u-) (89) 

where we have used total anitsymmetry of the structure constants. This 
completes our symmetrization procedure for ideal Yang-Mills fluids. 

In closing, we remark that the simple structure of the flux-freezing con- 
straints allows for explicit three-parametrizations 

for each Lie-component of the magnetic field, where the three-parametrized 
Lie-component kx satisfies k^^c = 0, which are such that cx = 0. By suitable 
choice of coordinates, = (1,0,0,0) so that ka = {0,ka)- Our foregoing 
analysis shows that in the 5 -|- 3N system parametrization V = {va, ka, /, T) 
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is such that Friedrichs' Properties CI and CII are satisfied (with the same 
Wa — {va,ha,f — TS,T) as before, but now in this reduced parametriza- 
tion and with respect to all variations in this parametrization) . Thus, in 
the case of such simple constraints, we could resort to pursuing Friedrichs' 
symmetrization procedure. We are then left with the task of proving that 
the symmetric hyperbolic system as defined by Friedrichs symmetrization 
procedure actually solves the entire system, now overdetermined by + 1. 
Further dependency relations must be invoked, namely those resulting from 
antisymmetry of hj^aUb]- We have not set out to pursue this, to see if such 
specific treatment results in shorter arguments, although this should lead to 
a successful argument as well. 
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